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ABSTRACT 

A  computer  program  was  written  for  the  HP  9845  and  the  HP  9816  to 
calculate  the  directivity  index  directly  from  the  beam  pattern.  A  graphics 
tablet  and  stylus  is  used  to  digitize  data  points  used  to  compute  the 
directivity  index.  Graphical  output  is  produced  on  a  printer  and  plotter. 
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INTRODUCTION 

This  program  was  developed  to  enable  the  user  to  digitally  input  beam 
pattern  data  points  into  a  computer  and  calculate  the  directivity  index  from 
this  data.  It  eliminates  the  slow  and  monotonous  task  of  visually  determining 
the  polar  coordinates  of  the  data  points  and  manually  keying  them  into  the 
computer  for  processing. 

Prior  to  the  creation  of  this  program,  two  methods  were  used  to  find  the 
directivity  index  of  a  beam  pattern.  The  first  method  involved  extracting 
data  from  the  beam  pattern  by  reading  a  data  point  at  specific,  regular 
intervals.  This  method  was  slow  and  less  accurate  due  to  the  fact  that  the 
operator  introduced  error  when  approximating  the  position  of  the  data  points. 
The  second  method  required  the  beam  pattern,  which  was  plotted  in  a  polar 
coordinate  system,  to  be  re-plotted  in  a  rectangular  coordinate  system.  This 
operation  was  performed  by  the  user.  Once  this  conversion  had  been  completed, 
a  mechanical  device  called  a  planimeter  was  used  to  calculate  the  directivity 
index.  Obviously,  this  method  was  extremely  slow.  Therefore,  this  program 
greatly  reduces  the  time  and  effort  required  to  calculate  the  directivity 
index  of  a  beam  pattern.  It  also  produces  a  more  accurate  result. 

The  program  requires  an  HP  9845  computer,  an  HP  9111A  graphics  tablet  and 
an  HP  9872A  plotter.  In  addition,  it  requires  that  the  Structured  Programming 
ROM  be  resident  in  the  HP  9845  computer.  Without  it,  certain  BASIC  statements 
cannot  be  executed. 

A  slightly  modified  version  of  this  program  is  available  for  the  HP  9816 
computer.  It  provides  the  same  capabilities  as  the  HP  9845  version  but 
utilizes  an  HP  7470  plotter  instead  of  an  HP  9872A  plotter.  The  HP  9816 
version  also  requires  that  the  BASIC  extensions  version  2.1  be  present. 
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DESCRIPTION  OF  PROGRAM 

The  "CALCDI"  program  is  responsible  for  converting  a  polar  plot 
representing  a  beam  pattern  data  into  a  finite  number  of  polar  coordinate  data 
points.  It  is  from  these  data  points  that  the  directivity  index  is  computed. 

The  program  is  divided  into  subprograms  which  perform  the  functions  of 
initializing  the  graphics  tablet,  digitizing  a  directivity  pattern, 
calculating  the  directivity  index,  and  producing  a  hard  copy  output. 
Subprograms  are  also  employed  to  perform  conversions  between  rectangular  and 
polar  coordinate  systems  and  perform  functions  not  found  in  HP  BASIC.  These 
programming  techniques  result  in  a  program  that  is  easy  to  understand  and 
modify.  Subprogram  modules  can  also  be  used  to  create  new  programs. 

PROBLEM  DEFINITION  AND  SOLUTION 

As  stated  in  the  introduction,  this  program's  purpose  is  to  allow  the 
operator  to  digitize  a  finite  number  of  data  points  from  a  directional 
response  pattern.  The  number  of  data  points  collected  is  dependent  upon  the 
incremental  angle  specified  by  the  operator.  After  collecting  the  data,  the 
beam  pattern's  directivity  index  calculation  is  performed  and  the  result  is 
displayed.  The  beam  pattern  is  plotted  on  the  CRT  to  allow  the  operator  to 
verify  that  no  errors  have  occurred  during  the  data  entry  process. 

Looking  at  the  program  from  an  overall  view,  it  can  be  divided  into  3 
distinct  steps.  The  first  step  involves  digitizing  a  beam  pattern  using  the 
graphics  tablet.  This  requires  converting  the  rectangular  coordinates  from 
the  graphics  tablet  to  the  polar  coordinates  required  by  the  directivity  index 
formula.  The  second  step  is  to  use  the  polar  coordinate  data  points  to 
compute  the  directivity  index.  This  entails  changing  the  formula  for 
calculating  the  directivity  index  into  a  procedure  able  to  be  executed  by  the 
computer.  The  directivity  index  formula,  an  integral,  is  translated  into  a 
summation.  Since  the  original  rectangular  coordinate  data  points  are  not 
saved,  the  polar  coordinate  data  points  must  be  converted  back  to  the 
rectangular  system  to  enable  them  to  be  displayed  on  the  CRT,  printer  or 
plotter. 


STEP  1 

DIGITIZING  A  DIRECTIONAL  RESPONSE  PATTERN 

A  beam,  or  directional  response,  pattern  is  digitized  by  placing  it  on 
the  surface  of  the  graphics  tablet  and  tracing  it  with  the  stylus.  Although 
the  beam  pattern  is  printed  on  polar  coordinate  paper,  X  and  Y  coordinates  are 
entered  into  the  computer  as  the  digitizing  process  occurs.  These  X  and  Y 
coordinates  must  be  converted  to  polar  form  for  use  by  the  directivity  index 
formula. 

Polar  coordinate  points  are  in  the  form  (Rho,  Theta)  with  Rho  equalling 
the  distance  from  the  origin  and  Theta  equalling  an  angle  in  the  range  0-360 
degrees.  Rho  and  Theta  must  be  derived  from  the  X  and  Y  coordinates.  This 
conversion  process  is  accomplished  by  the  following  subprograms: 
"Displacement,"  "Calculateradius,"  "Calcarctangent,"  "Determinequad,"  and 
Convertto360."  Each  of  these  subprograms  is  well  commented.  The  procedure 
for  the  conversion  is  also  presented  in  algorithmic  form. 
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STEP  2 

CALCULATE  THE  DIRECTIVITY  INDEX 

The  directivity  index  calculation  is  the  most  important  segment  of  the 
program.  From  the  polar  coordinate  data,  we  are  able  to  determine  the 
directivity  index.  This  computation  is  performed  by  the  subprogram  called 
"Calcdirindex." 

The  directivity  index  is  equal  to  ten  times  the  common  logarithm  of  the 
directivity  factor.  The  directivity  factor  is  defined  as  the  ratio  of  the 
intensity  in  a  reference  direction,  usually  the  axis,  to  the  intensity  average 
over  all  directions. 

The  directivity  factor  is  calculated  using  the  formula  as  presented  by 
Bobber  in  Underwater  Electroacoustic  teasurements  (Washington  D.C.,  NRL,  1970, 
P-  84) 

From  Bobber, 


directivity  factor 


n 

A9  .£ 

i=0 


Ai 


(1) 


where  a  e  equals  the  degrees  resolution  per  angular  interval,  n  equals  the 
number  of  sampling  points  per  180  degrees,  and 


.  sin  oi 


(2) 


In  equation  2,  en-  is  the  angular  measure  in  degrees  of  the  polar  coordinate 
data  point.  Also, 


l0-(db  down/10) 


(3) 


where  db  down  equals  the  difference  of  the  maximum  response  axis  and  the 
radius  of  the  polar  coordinate  data  point.  The  formula  for  sampling  the 
directivity  pattern  through  360  degrees  becomes. 


directivity  factor  =  - 

a  e 


(4) 
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Consider  the  case  of  a  e  =  5  degrees,  the  default  value  for  the  resolution  in 
''Calcdirindex."  Converting  a  e  to  radians. 


directivity  factor  = 


(5) 


Simpl ifying. 


directivity  factor 


229.2 


(6) 


The  directivity  factor  formula  of  Equation  6  is  transformed  into  a 
sequence  of  BASIC  statments.  "Calcdirindex"  simulates  the  summation  process 
by  utilizing  an  iterative  loop.  The  variable  "resolution"  corresponds  to 
a  e  and  229.2  is  a  numeric  constant. 

If  a  directivity  pattern  has  symmetry  about  the  equator  of  the  measuring 
sphere,  then  it  is  only  necessary  to  use  one  half  of  the  pattern  in  the 
computation.  However,  some  patterns  are  nearly  rotational ly  symmetric,  but 
not  exactly.  For  these  "near  symmetric  patterns,"  it  is  desirable  to  average 
the  two  directivity  factors  as  calculated  for  the  two  halves  of  the  pattern. 
Therefore,  it  must  be  noted  that  this  program  assumes  all  beam  patterns  are 
nearly  rotational ly  symmetric  and  thus  performs  the  summation  over  the  full 
360  degrees. 


STEP  3 

DISPLAY  THE  DIRECTIVITY  PATTERN  ON  THE  OUTPUT  DEVICE 

The  polar  coordinate  data  is  converted  to  X  and  Y  coordinates  so  that  it 
may  be  output  to  the  CRT,  printer,  and  plotter.  This  is  accomplished  by  the 
subprograms  called  "Displaypattern,"  "Printcopy,"  and  "Plotcopy." 
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Pseudo-code 
Mai n  Program 


1.  Clear  the  CRT 

a.  Clear  the  alphanumeric  raster 

b.  Clear  the  graphics  raster 

2.  Print  Identifying  message 

3.  Wait  for  person  to  read  message 

a.  Continue  when  operator  hits  "Continue"  key 

4.  Digitize  the  directivity  pattern  using  the  9111A  Graphics  Tablet 

5.  Calculate  the  directivity  index  from  the  polar  coordinate  values  entered 

from  the  graphics  tablet 

6.  Draw  the  directivity  pattern  on  the  CRT  and  print  the  directivity  index 

7.  Output  directivity  pattern  to  printer  with  option  of  suppressing  output 

8.  Output  directivity  pattern  to  plotter  with  option  of  suppressing  output 

9.  Clear  the  CRT 

a.  Clear  the  alphanumeric  raster 

b.  Clear  graphics  raster 

10.  See  if  operator  has  more  directivity  patterns  to  digitize 

If  choice  =  yes  then 

Repeat  procedure  from  step  4 

Else 

End  of  program 
End  if 
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Pseudo  Code 

Algorithm  to  Convert  from  the  Rectangular 
to  the  Polar  Coordinate  System 


1.  Input  the  X  and  Y  coordinates 

2.  Determine  its  X  and  Y  displacements  from  the  origin  in  user  units 

3.  Calculate  Rho  (Radius)  using  the  X  and  Y  displacements  with  Rho  being 

measured  in  user  units 

4.  Compute  the  principle  angle  given  the  X  and  Y  displacements  from  the  origin 

5.  Determine  which  quadrant  the  point  lies  in  based  upon  its  relative  position 

to  the  origin 

6.  Given  the  principle  angle  and  the  quadrant,  produce  an  angle  in  the  range 

0-360  degrees 
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VARIABLES 

originx  -  the  x-coordinate  of  the  origin  of  the  directivity  pattern  in 
digitizing  units. 

originy  -  the  y-coordinate  of  the  origin  of  the  directivity  pattern  in 
digitizing  units. 

zerodegfsx  -  the  x-coordinate  of  the  point  representing  zero  degrees  full 

scale  on  the  directivity  pattern.  It  is  measured  in  digitizing 
units. 

zerodegfsy  -  the  y-coordinate  of  the  point  representing  zero  degrees  full 

scale  on  the  directivity  pattern.  It  is  measured  in  digitizing 
units. 

mrax  -  the  x-coordinate  of  the  point  representing  the  maximum  response  axis. 

It  is  measured  in  digitizing  units. 

mray  -  the  y-coordinate  of  the  point  representing  the  maximum  response  axis. 

It  is  measured  in  digitizing  units. 

gtabuperdb  -  number  of  digitizing  units  (graphic  tablet  units)  per  decibel. 
This  is  used  as  a  scaling  or  conversion  factor. 

dbfullscale  -  number  of  decibels  from  the  origin  to  zero  degrees  full  scale. 

radiusvalues  (*)  -  an  array  representing  the  radii  of  the  polar  coordinate 

data  points. 

angle  (*)  -  an  array  representing  the  angles  of  the  polar  coordinate  data 
points. 

mra  -  the  maximum  response  axis  value. 

resolution  -  the  incremental  angle,  which  determines  the  number  of  data  points 
gathered. 

dirindex  -  the  value  of  the  directivity  index. 

skew  -  an  offset  angle  to  compensate  for  directivity  patterns  positioned 
crookedly  on  the  graphics  tablet  digitizing  area. 
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*********************************************************************** 


*  * 

*  AUTOMATIC  DIRECTIVITY  INDEX  CALCULATION  PROGRAM  * 

*  by  JEHU  ROOKARD  * 

*  * 

*  Created  2  AUG  1984  * 

*  * 


*********************************************************************** 

THIS  PROGRAM  DIGITIZES  A  DIRECTIVITY  PATTERN  AND  CALCULATES  THE 
DIRECTIVITY  INDEX 


Main:  !  MAIN  DIRECTIVITY  PATTERN  PROGRAM 

OPTION  BASE  0 

COM  Originx,Originy,Zerodegfsx,Zerodegfsy,Mrax,Mray 

COM  Gtabuperdb, Dbful  1  seal  e 

COM  Radi  usual ues<360> , Angl  e<360) 

COM  Mra 
COM  Resolution 
COM  Dirindex 
COM  Skew 

CALL  Clearscreen 

PLOTTER  IS  13, "GRAPHICS"  !  ACTIVATE  CRT  GRAPHICS  RASTER 
GRAPHICS 

GCLEAR  !  CLEAR  THE  GRAPHICS  RASTER 

PLOTTER  13  IS  OFF  !  DEACTIVATE  CRT  GRAPHICS  RASTER 

PRINT  "THIS  PROGRAM  ENTERS  A  DIRECTIVITY  PATTERN  FROM  THE  GRAPHICS" 

PRINT  "TABLET  AND  CALCULATES  THE  DIRECTIVITY  INDEX" 

INPUT  "PRESS  " CONT '  WHEN  READY  TO  PROCEED" , Cont  $ 

CALL  Enterapattern  !  ENTER  A  DIRECTIVITY  PATTERN  USING  THE  9111A 

!  GRAPHICS  TABLET 

CALL  Calcdirindex  !  CALCULATE  THE  DIRECTIVITY  INDEX 

CALL  Di spl aypattern  !  DRAW  THE  DIRECTIVITY  PATTERN  ON  THE  CRT 

CALL  Printcopy  !  MAKE  HARD  COPY  ON  PRINTER 

CALL  Plotcopy  !  MAKE  HARD  COPY  ON  PLOTTER 

PLOTTER  IS  13, "GRAPHICS"  !  CLEAR  THE  GRAPHICS  RASTER 

GCLEAR 

EXIT  GRAPHICS 

PLOTTER  13  IS  OFF  !  DEACTIVATE  CRT  GRAPHICS  RASTER 

CALL  Clearscreen 

Yesno$="Y" 

INPUT  "ENTER  ANOTHER  PATTERN  -  Y/N.  DEFAULT  IS  ' Y ' " , Yesno* 

IF  Yesno$="Y"  THEN 

PLOTTER  IS  13, "GRAPHICS"  !  CLEAR  THE  GRAPHICS  RASTER 
GCLEAR 

PLOTTER  13  IS  OFF  -  !  DEACTIVATE  THE  CRT  GRAPHICS  RASTER 

GOTO  310 

ELSE 

END 


Enterapattern:  !  PROCEDURE  TO  ENTER  A  DIRECTIVITY  PATTERN  FROM  THE  9111A 

!  GRAPHICS  TABLET 
SUB  Ent erapat t ern 

CALL  Clearscreen  !  CLEAR  THE  ALPHA  RASTER 

PLOTTER  IS  7,6, "9372A"  !  THE  9111A  GRAPHICS  TABLET  CAN  ONLY  BE  ACCESSED 

!  THROUGH  THE  GRAPHICS  ROM.  ALTHOUGH  IT  IS  A 
!  GRAPHICS  TABLET,  THE  COMPUTER  TREATS  IT  AS  A 
!  PLOTTER 

OUTPUT  706; "IN"  !  INITIALIZE  THE  9111A  GRAPHICS  TABLET 

CALL  Setupaxes  !  GET  THE  NECESSARY  POINTS  TO  PROPERLY  SCALE 

!  THE  DIGITIZING  AREA  AND  ALL  POINTS  ENTERED 
!  FROM  IT 


660  CALL  Entercurve  !  DIGITIZE  THE  ACTUAL  DIRECTIVITY  PATTERN 

670  SUBEND 
680  ! 

690  ! 

700  ! 

710  Set upaxes :  !  ENTER  THE  NECESSARY  SCALING  VALUES  FROM  THE  GRAPHICS 

720  !  TABLET  SO  THAT  A  BASIS  FOR  ENTERING  ALL  SUCCEEDING 

730  !  POINTS  IS  ESTABLISHED 

740  !  A  METHOD  (SCALE)  FOR  CONVERTING  BETWEEN  DIGITIZING  UNITS 

750  !  AND  USER  UNITS  IS  SET  UP 

760  SUB  Setupaxes 

770  OPTION  BASE  0 

780  COM  Or i gi nx, Ori gi ny , Zerodegfsx, Zerodegfsy, Mrax, Mray 

790  COM  Gtabuperdb,  Dbful  1  seal  e 

300  COM  Radi  usual ues(360) , Angl e(360) 

810  COM  Mr a 

820  COM  Resolution 

330  COM  Dirindex 

840  COM  Skew 

850  Dbful 1  seal e  =  50 

860  INPUT  "ENTER  NO.  OF  db'S  FROM  ORIGIN  TO  ZERO  DEGREES  FULL  SCALE.  DEFAULT  I 

S  50" , Dbful 1  seal e 

870  IF  Dbful 1  seal e>5@  THEN 

380  PRINT  "NO  VALUES  GREATER  THAN  50  ALLOWED.  PLEASE  RE-ENTER" 

890  GOTO  850 
900  ELSE 
910  END  IF 
920  ! 

930  Resolutions 

940  INPUT  "ENTER  INCREMENTAL  ANGLE.  DEFAULT  IS  5  DEGREES ", Reso 1 ut i on 
950  Resol ut i on=INT(Resol ut i on) 

960  IF  -Reso  1  ut  i  on<  1  THEN 

970  PRINT  "INCREMENTAL  ANGLE  MUST  BE  GREATER  THAN  OR  EQUAL  TO  1" 

980  GOTO  930 
990  ELSE 
1000  END  IF 
1010  ! 

1020  DISP  "DIGITIZE  ORIGIN" 

1030  OUTPUT  706; "SG"  !  SET  TO  SINGLE  POINT  MODE 

1040  CALL  Status  !  WAIT  UNTIL  GRAPHICS  TABLET  IS  READY 

1050  OUTPUT  706; "OD"  !  INSTRUCT  GRAPHICS  TABLET  TO  TRANSMIT  POINT 

1060  ENTER  706; Or i gi nx, Or i gi ny  !  GET  THE  POINT  FROM  THE  GRAPHICS  TABLET 

1070  OUTPUT  706; "BP24, 125, 4"  !  INTRUCT  GRAPHICS  TABLET  TO  BEEP 

1080  ! 

1090  DISP  “DIGITIZE  ZERO  DEGREE  FULL  SCALE" 

1100  CALL  Status 
1110  OUTPUT  706; "OD" 

1120  ENTER  706 ; Zerodegfsx , Zerodegfsy 
1130  OUTPUT  706; "BP" 

1140  ! 

1150  DISP  "DIGITIZE  MAXIMUM  RESPONSE  AXIS" 

1160  CALL  Status 
1170  OUTPUT  706; "OD" 

1180  ENTER  706; Mrax, Mray 
1190  OUTPUT  706; "BP" 

1200  ! 

1210  !  CALCULATE  THE  NUMBER  OF  DIGITIZING  UNITS  PER  DECIBEL  BY  DIVIDING 

1220  !  THE  RADIUS  (Rho)  AT  FULL  SCALE  BY  THE  NUMBER  OF  DECIBELS  FULL  SCALE 

1230  Differencex=ABS(Zerodegfsx-Originx) 

1240  Di f fere ncey=ABS(Zerodegfsy- Or i gi ny) 

1250  CALL  Calculateradius(Differencex,Differencey, Rho_zerodegf s ) 

1260  Gt  abuperdb  =  Rho_zerodegf sxDbf u 1 1  sc  al e 
1270  ! 

1280  !  CALCULATE  THE  MAXIMUM  RESPONSE  AXIS  VALUE  BY  DIVIDING  THE  VALUE  FOR 

1290  !  THE  MAXIMUM  RESPONSE  AXIS  IN  DIGITIZING  UNITS  BY  THE  NUMBER  OF 

1300  !  DIGITIZING  UNITS  PER  DECIBEL 


1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660. 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 

1820 

1330 

1840 

1350 

I860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

I960 


!  <  ie.  CONVERT  FROM  DIGITIZING  UNITS  TO  USER  UNITS  > 
Differencex=ABS(Mrax-Originx) 

Di fferencey=ABS(Mray-Or i gi  ny) 

CALL  Calculateradius(Differencex,Differencey,Mra) 

Mra=MraxGt  abuperdb 

!  DETERMINE  THE  "SKEW  ANGLE"  TO  CORRECT  FOR  ERRORS  THAT  WILL  RESULT  FROM 
!  THE  ORIGINAL  DIRECTIVITY  PATTERN  BEING  POSITIONED  CROOKED  OR  SKEWED  ON 
!  THE  GRAPHICS  TABLET.  THIS  CORRECTION  FACTOR  ALLOWS  THE  USER  TO  PLACE  THE 
!  ORIGINAL  DIRECTIVITY  PATTERN  IN  any  POSITION  inside  THE  GRAPHICS  TABLET 
!  DIGITIZING  AREA  AND  STILL  GATHER  ACCURATE  DATA.  IT  ALSO  RESULTS  IN  A 
!  TRUE  POLAR  CO-ORDINATE  SYSTEM 

CALL  DispIacement(Zerodegfsx,Zerodegfsy,Xualue,Y'v'alue) 

CALL  Calc arctangent<Xoalue,Yualue,Rel at ioe angle) 

CALL  Determi  nequad ( Zerodegf sx , Zerodegf sy, Quadrant  ) 

CALL  Co noertto360(Rel at i  oe angle, Quad rant, Absolute angle) 

Skew=Absolute angle 

PLOTTER  7,6  IS  OFF  !  DEACTIVATE  GRAPHICS  TABLET 

!  <  ie.  DON'T  SEND  GRAPHICS  COMMANDS  TO  IT  ) 

SUBEND 


Entercurue.'  !  ENTER  A  DIRECTIVITY  PATTERN  USING  THE  9 1 1 1 A  GRAPHICS  TABLET 
!  THE  USER  MUST  START  DIGITIZING  FROM  ZERO  DEFREES  AND 
!  CONTINUE  UNTIL  HE  REACHES  360  DEGREES.  'Resolution' 

!  DETERMINES  THE  NUMBER  OF  DATA  POINTS  ACCEPTED. 

!  360x'Resolut ion'  SECTORS  ARE  CREATED  AND  A  DATA  POINT 
!  IS  ENTERED  FOR  EACH  SECTOR  AND  STORED  IN  ITS  RESPECTIVE 
!  ARRAY  ELEMENT 
SUB  Entercuroe 
OPTION  BASE  0 

COM  Originx,Originy,Zerodegfsx,Zerodegfsy,Mrax,Mray 

COM  Gtabuperdb, Dbful 1  seal  e 

COM  Radiusualues<360),Angle<360) 

COM  Mra 
COM  Reso 1 ut i on 
COM  Dirindex 
COM  Skew 

DEG  !  DEGREES  MODE  < NOT  RADIANS) 

! 

!  "ENTERCURVE"  EXPECTS  TO  HAVE  A  "COMMUNICATIONS  CHANNEL"  WITH  THE 
!  GRAPHICS  TABLET  ALREADY  ESTABLISHED.  IT  ALSO  EXPECTS  THE  GRAPHICS 
!  TABLET  TO  HAVE  BEEN  INITIALIZED.  <  AS  IN  "SETUPAXES"  ) 

!  "ENTERCURVE  WILL  RE-ACTIVATE  THE  GRAPHICS  TABLET  IF  IT  HAS  PREVIOUSLY 
!  BEEN  DEACTIVATED 

j 

PLOTTER  7,6  IS  ON 

OUTPUT  706; "CN"  !  SET  GRAPHICS  TABLET  TO  CONTINUOUS 

!  SAMPLING  MODE 

DISP  "START  DIGITIZING  FROM  ZERO  DEGREES  AND  PROCEED  COUNTER-CLOCKWISE" 

FOR  Angl ecounter-0  TO  360-Resolution  STEP  Resolution 

CALL  Status  !  WAIT  UNTIL  GRAPHICS  TABLET  IS  READY 

OUTPUT  706; "OD" 

ENTER  706; Beamf ormx , Beamf ormy  !  GET  DATA  POINT 

! 

!  CONVERT  THE  X  AND  Y  COORDINATES  TO  POLAR  COORDINATES  BY  CALLING  THE 
!  FOLLOWING  SUBPROGRAMS  IN  THE  PROPER  SEQUENCE 

! 

CALL  Displacement<Beamformx,Beamformy,Xdisp1acement,Ydisplacement) 

CALL  CalculateradiusCXdispl acement  ,Ydispl acement, Radius) 

CALL  Calc ar ct an gent<Xdisplacement,Ydisplacement,Rel at  ioe angle) 

CALL  Determinequad(Beamformx,Beamformy, Quadrant ) 

CALL  Conoertto360(Relati oeangl e, Quadrant ,  Absol  uteangl  e) 

! 

!  TEST  TO  SEE  IF  THE  POINT  LIES  WITHIN  THE  SECTOR.  IF  IT  DOES,  THEN  STORE 


1970  !  IT  IN  ITS  RESPECTIVE  ARRAY  ELEMENT.  REMEMBER  TO  CORRECT  FOR  ANY  "SKEW" 

1980  !  ERRORS  THAT  MAY  EXIST 

1990  ! 

2000  IF  < Abso 1 ut eang 1 e-Sk eu >Ang 1 ec ount er )  AND  < Abso 1 ut eang 1 e-Sk ew< Ang  1  ec ount er  + 
Resol ut i on )  THEN 
2010  ELSE 
2020  ! 

2030  !  IF  'Skew'  IS  NEGATIVE  THEN  PERFORM  THE  TEST  BUT  WITH  360  DEGREES  ADDED 

2040  !  TO  MAKE  THE  ANGLE  POSITIVE 

2050  ! 

2060  IF  (Absol uteangl e+360-Skew>Angl ecounter)  AND  < Abso 1 ut eang 1 e+360-Sk ew< Ang 1 e 
counter+Resol ut i on)  THEN 


2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2288 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 


Absoluteang1e=Abso1uteangle+360 

ELSE 

GOTO  1830 
END  IF 
END  IF 

DISP  "ANGLE  =" ; Absol uteangl e-Skew, "db  DOWN  = " ; Mra-Rad i us 
Radi  usual ues<Angl ecounter)=R ad i us 
Angle<Anglecounter)=Absoluteangle-Skew 
BEEP 

NEXT  Anglecounter 


PLOTTER  7,6  IS  OFF  !  DEACTIVATE  GRAPHICS  TABLET 

SUBEND 


D 


sp  1 ac  ement : 


GIVEN  A  POINT'S  X  AND  Y  CO-ORDINATES,  DETERMINE  IT'S 
RELATIVE  X  AND  Y  DISPLACEMENTS  FROM  THE  ORIGIN 
' Xval ue '  IS  THE  X-CO-ORD I  NATE  OF  THE  POINT  IN  GRAPHIC 
TABLET  UNITS- 

'Yvalue'  IS  THE  Y-CO-ORDINATE  OF  THE  POINT  IN  GRAPHIC 
TABLET  UNITS 

' Xd i sp 1 ac  ement '  IS  THE  RELATIVE  X-D I SPLACEMENT  FROM  THE 
ORIGIN  IN  USER  DEFINED  UNITS 

' Ydi spl acement '  IS  THE  RELATIVE  Y-D I SPLACEMENT  FROM  THE 
ORIGIN  IN  USER  DEFINED  UNITS 


SUE  Displacement<Xvalue,Yvalue,Xdisp1ac ement, Ydisplac ement) 
COM  Originx,Originy,Zerodegfsx,Zerodegfsy,Mrax,Mray 
COM  Gtabuperdb, Dbful 1  seal  e 

Xdisplacement=<Xva1ue-0riginx)xGtabuperdb 

Ydisplacement=<Yva1ue-0riginy)xGtabuperdb 

SUBEND 


Cal cul ateradi us:  !  CALCULATE  THE  RADIUS  USING  THE  PYTHAGAREAN  THEOREM 

!  <  D I  STANCE  FORMULA) 

SUE  Cal cul ateradi us <X, Y, R) 

R=s$QR(Xa2+Ya2) 

SUBEND 


Cal  care  tangent:  !  CALCULATE  THE  VALUE  OF  THE  ARCTANGENT  GIVEN  'X'  AND  'Y' 

!  ALSO  RESOLVE  ANY  INVALID  INPUT  PARAMETERS  TO  THE 
!  ARCTANGENT  FUNCTION  BY  SUBSTITUTING  THE  CORRECT  VALUE 
SUB  Cal carctangent <X, Y, Angl e) 

DEG  !  DEGREES  MODE  < NOT  RADIANS) 

ON  ERROR  GOTO  Recovery 
Angl e=ATN<Y/X) 

SUBEXIT 

Recovery:  !  ROUTINE  TO  RECOVER  INVALID  INPUT  PARAMETERS  TO 

!  ARCTANGENT  FUNCTION.  <ie.  DENOMINATOR  CAN'T  EQUAL 
!  ZERO) 


2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3100 

3110 

3120 

3130 

3140 

3150 

3160 

3170 

3180 

3190 

3200 

3210 

3220 

3230 

3240 

3250 

3260 


!  ANGLE=ATN< Y/X) 
IF  ERRN=3 1  THEN 


SELECT  SGN< Y) 


CRSE  1 


RSSUMES  X  DOES  NOT  0 
IF  DIVISION  BY  ZERO  ERROR  THEN  CORRECT  THE  ERROR 
BY  ASSIGNING  THE  CORRECT  VALUE  OR  SOME  ARBITRARY 
VALUE  TO  ' ANGLE'' 


IF 

(  i  « 


A  NUMBER  IS  EITHER  +  ,-,  OR  ZERO,  SO  DETERMINE  THE 
CORRECT  VALUE  FOR  THE  3  CASES  OF  '  Y' 

'Y'  IS  POSITIVE  THEN  THE  CORRECT  ANSWER  IS  90  DEGREES 
.  THE  POINT  LIES  SOMEWHERE  ON  THE  Y-AXIS,  BUT  ABOVE 


Ang  1 e=0 
CASE  0 


Ang 1 e  =  45 1 


CASE  -1 


Ang 1 e=0 

END  SELECT 

SUBEXIT 

ELSE 


THE  X-AXIS) 

0  PRODUCES  A  VALUE  OF  90  DEGREES  AFTER  BEING  PROCESSED 
BY  " CONVERTTO3601* 

IF  'Y'  IS  ZERO  THEN  THE  POINT  IS  AT  THE  ORIGIN,  SO 
ARBITRARILY  ASSIGN  AN  ANGLE  OF  ZERO  DEGREES.  THE  VALUE 
OF  THE  RHO  (RADIUS)  WILL  BE  CALCULATED  AS  ZERO  AND  IT  WILL 
BE  OBVIOUS  TO  THE  USER  THAT  THE  POINT  WAS  AT  THE  ORIGIN 
451  PRODUCES  A  VERY  LARGE  VALUE  AFTER  BEING  PROCESSED  BY 
"CONVERTTO360".  THIS  RESULTS  IN  " ENTERCURVE "  REJECTING 
ANY  POINTS  THAT  LIE  AT  THE  ORIGIN  BECAUSE  THEY  WILL  BE 
OUT  OF  RANGE 

IF  'Y'  IS  NEGATIVE  THEN  THE  CORRECT  ANSWER  IS  270  DEGREES 
<ie.  THE  POINT  LIES  SOMEWHERE  ON  THE  Y-AXIS,  BUT  BELOW  THE 
X-AXIS) 

0  PRODUCES  A  VALUE  OF  270  DEGREES  AFTER  BEING  PROCESSED 
BY  " CONVERTTO360" 

NO  NEED  FOR  A  'CASE  ELSE'  STATEMENT  AS  A  NUMBER  IF  EITHER 
+,-  OR  ZERO 

END  OF  ERROR  HANDLER  ROUTINE 

IF  THIS  ISN'T  A  DIVISION  BY  ZERO  ERROR,  THEN  DISPLAY  THE 
ERROR  MESSAGE  AND  PAUSE  PROGRAM  EXECUTION 


DISP  ERRM$ 
PAUSE 
END  IF 
SUBEND 


Determi nequad:  !  DETERMINE  WHICH  QUADRANT  THE  POINT  LIES  IN  BASED  UPON 

!  WHETHER  THE  POINT  LIES  TO  THE  LEFT  OR  RIGHT  OF  THE 
!  Y-AXIS  AND  ABOVE  OR  BELOW  THE  X-AXIS 


SUB  Deter minequ ad <Xva1ue,Yvalue, Quadrant) 

COM  Originx,Originy,Zerodegfsx,Zerodegfsy,Mrax,Mray 


IF  Xval ue>Or i gi nx  THEN 
IF  Yval ue>=Or i gi ny  THEN 
Quadrant=l 
ELSE 


!  IN  QUADRANT  1  OR  4  THEN 
!  IN  QUADRANT  1 


Quadrant  =4 
END  IF 
ELSE 

IF  Yval ue >Or i g i ny  THEN 

Quadrant  =2 

ELSE 

IF  Xval ue=Or i gi nx  THEN 

Quadrant  =4 

ELSE 


!  IN  QUADRANT  2  OR  3  OR  ON  Y-AXIS  THEN 
!  IN  QUADRANT  2 

! IN  QUADRANT  4  (  ON  Y-AXIS  ) 


Quadrant  =3 


END  IF 
END  IF 
END  IF 
SUBEND 


Convert- 1 o360 :  !  GIVEN  THE  PRINCIPAL  ANGLE  AND  THE  QUADRANT  PRODUCE 

!  AN  ANGULAR  VALUE  IN  THE  RANGE  0  TO  360  DEGREES 
!  'Rel at i veangl e'  IS  THE  PRINCIPAL  ANGLE 
!  'Quadrant'  IS  IN  THE  RANGE  1-4  <  FROM  CARTESIAN 


3270  !  CO-ORDINATE  SYSTEM  > 

32S0  !  "flbsol uteangl e'  IS  THE  RESULTANT  ANGULAR  VALUE  IN 

3290  '  !  THE  RANGE  0-360  DEGREES 

3300  SUE  Convertto360<Rel  at  i  veangl  e,  Quadrant ,  Absol  uteangl  e) 

3310  SELECT  Quadrant 

3320  CASE  1  !  FIRST  QUADRANT  0  <=  Re  1  at i veang 1 e  <  90 

3330  Absolute angle=Rel at ive angle +270 

3340  CASE  2  !  SECOND  QUADRANT  90  <=  Re  1  at i veang 1 e  <  180 

3350  Absolute an gle=Re1 at ive an gle+ 90 

3360  CASE  3  !  THIRD  QUADRANT  180  <=  Re  1  at i veang 1 e  <  270 

3370  Abso 1 ut  eang 1 e  =  Re 1  at i veang 1 e+90 

3380  CASE  4  !  FOURTH  QUADRANT  270  <=  Re  1  at i ueang 1 e  <  360 

3390  Abso 1 ut  eang 1 e=Re lativeangle+270 

3400  END  SELECT 

3410  SUBEND 

3420  ! 

3430  ! 

3440  ! 

3450  Clearscreen:  !  CLEAR  THE  ALPHA  RASTER  BY  PRINTING  24  BLANK  LINES 

3460  SUB  Clearscreen 

3470  FOR  L i ne=  1  TO  24 

3480  PRINT 

3490  NEXT  Line 

3500  SUBEND 

3510  ! 

3520  ! 

3530  ! 

3540  Di spl aypattern:  !  DRAW  THE  DIRECTIVITY  PATTERN  ON  THE  GRAPHICS  RASTER 

3550  !  AND  LABEL  THE  DIRECTIVITY  INDEX 

3560  SUB  D i sp 1 aypat t ern  . 

3570  OPTION  BASE  0 

3580  COM  Or i gi nx, Or i gi ny, Zerodegfsx, Zerodegfsy , Mrax, Mray 

3590  COM  Gt abuperdb , Dbf u 1  1  sc al  e 

3600  COM  Radi  usual ues<360) , Angl e<360) 

3610  COM  Mra 
3620  COM  Resolution 
3630  COM  Dirindex 

3640  DEG  !  DEGREES  MODE  <  NOT  RADIANS  > 

3650  PLOTTER  IS  1 3 , " GRAPH  I CS "  !  ACTIVATE  THE  CRT  GRAPHICS  RASTER 

3660  GRAPHICS  !  ENABLE  GRAPHICS 

3670  GCLEAR  !  CLEAR  GRAPHICS  RASTER 

3690  LOCATE  20,95,0,100 

3690  MOVE  37,2 

3700  CSI2E  3.2 

3710  LABEL  "Directivity  index  =" ; DROUNDCDi r i ndex,  6) 

3720  SCALE  -Dbful 1  sc al e* . 75 , Dbf u 1 1  sc al e*. 75, -Dbful 1  seal e, Dbf ul 1  seal  e 
3730  FRAME 

3740  AXES  10,10,0,0,1,1 

3750  FOR  Ang 1 ec ount er=@  TO  360-Resolution  STEP  Resolution 
3760  X=Radi  usual  ues  ( Ang  1  ec  ount  er  >*COS<  Ang  1  e  <  Ang  1  ecount  er  )  +90  !> 

3770  Y=Radi usual ues ( Ang 1 ecounter >*SIN(Angl e<Anglecounter>+90> 

3780  IF  Anglecounter=0  THEN 

3790  MOVE  X , Y 

3800  ELSE 

3810  DRAW  X , Y 

3820  END  IF 

3830  NEXT  Anglecounter 

3840  !  COMPLETE  THE  DIRECTIVITY  PATTERN  PLOT  BY  DRAWING  A  LINE  FROM  THE 

3850  !  LAST  POINT  TO  THE  FIRST  POINT. 

3860  DRAW  Radi  usual ues<0)*COS<Angle<0>+90), Rad i usual ues<0>*SIN<Angle<0)+90) 
3870  PLOTTER  13  IS  OFF  !  DEACTIVATE  THE  CRT  GRAPHICS  RASTER 

3380  SUBEND 
3890  ! 

3900  ! 

3910  ! 

3920  Cal edi r i ndex:  !  CALCULATE  THE  DIRECTIVITY  INDEX  GIVEN  THE  RADIUS  <Rho> 


3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

0) 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 

4220 

4230 

4240 

4250 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4430 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 


!  AND  THE  ANGLE  <  0-360  DEGREES  ) 

SUB  Calcdirindex 
OPTION  BASE  0 

COM  Ori gi nx, Ori gi ny, Zerodegfsx, Zerodegfsy, Mrax, Mray 

COM  Gtabuperdb, Dbf ul 1  seal e 

COM  Radiusualues(360),Ang!e<360) 

COM  Mra 

COM  Resolution 

COM  Di ri ndex 

DEG  !  DEGREES  MODE,  NOT  RADIANS 

Sum  =  0 

Ang 1 ec  ount  er  =  0 

FOR  Angl ecounter=0  TO  360-Resolution  STEP  Resolution 

Temp=ABS<SIN<Angle(Anglecounter)))*10/'<-<Mra-Radiusua1ues<Anglecounter))xl 

Sum  =  T  emp  +  Sum 

NEXT  Anglecounten 

Di rf actor=229.2/<Resolution*Sum) 

Dirindex=10*LGT(Dirfactor) 

SUBEND 


Pr i ntcopy :  !  DUMP  THE  CONTENTS  OF  THE  GRAPHICS  RASTER  TO  THE  PRINTER 

!  ONLY  THE  DIRECTIVITY  PATTERN  AND  ITS  RESPECTIVE  LABELS  WILL 
!  APPEAR  ON  THE  PRINTER  OUTPUT 
SUB  Printcopy 
Yesno*="N" 


INPUT  "DUMP  DIRECTIVITY  PATTERN  TO  PRINTER  -  YxN.  DEFAULT  IS  ■'N'",Yesno$ 
IF  Yesnol="Y"  THEN 
PLOTTER  13  IS  ON 
PRINTER  IS  0 
DUMP  GRAPHICS 
PRINT  CHR# < 12 ) 

PRINTER  IS  16 
PLOTTER  13  IS  OFF 
ELSE 
END  IF 
SUBEND 


!  OUTPUT  DEVICE  =  PRINTER 

!  ADVANCE  TO  TOP  OF  NEXT  SHEET  OF  PAPER 
!  OUTPUT  DEVICE  =  CRT 


Status:  !  INTERROGATE  THE  9111A  GRAPHICS  TABLET  TO  SEE  IF  IT'S  READY 

!  READ  THE  STATUS  WORD  AND  TEST  BIT  2.  THE  GRAPHICS  TABLET  IS 
!  READY  WHEN  BIT  2  =  1 
SUB  Status 
STATUS  7,6;S 

IF  BIT<S,2>=0  THEN  4380  - 
SUBEND 


Plotcopy:  !  PLOT  THE  DIRECTIVITY  PATTERN  ON  THE  PLOTTER  AND  LABEL  THE 

!  DIRECTIVITY  INDEX 
SUB  Plotcopy 
OPTION  BASE  0 

COM  Or  i  ginx,  Ori  giny,  Zerodegfsx,  Zerodegfsy,  Mr  ax,  Mr  ay 

COM  Gt abuperdb , Dbf u 1 1  sc al e 

COM  Radi  usual ues< 360) , Angl e(360) 

COM  Mra 

COM  Resolution 

COM  Dirindex 

DEG  !  DEGREES  MODE  <NOT  RADIANS) 

Yesno*=“N" 

D i mf 1 ag=0 

INPUT  "DUMP  DIRECTIVITY  PATTERN  TO  PLOTTER  -  Y/N.  DEFAULT  IS  'N/",Yesnof 


4580  IF  Yesno$=" Y"  THEN 
4590  IF  D i mf 1 ag= 1  THEN  4630 
4600  DIM  T emp$ <  2  ) C  35  ] 

4610  Dimflag=l 

4620  PLOTTER  IS  7,5,"9872A“ 

4630  Temp$<0>="  " 

4640  Temp$a>  =  "  " 

4650  Temp$<2)="  " 

4660  INPUT  "ENTER  LABEL  FOR  FIRST  LINE" , Temp$<0) 

4670  INPUT  "ENTER  LABEL  FOR  SECOND  LINE", T emp$ <  1  ) 

4680  INPUT  "ENTER  DATE  (OPTIONAL Temp*<2) 

4690  LOCATE  0 , R AT  I  0* 1 00 , 0 , 1 00 

4700  FRAME 

4710  LORG  2 

4720  LDIR  90 

4730  CSIZE  2.6 

4740  MOVE  120,4  !  DRAW  LABEL  1 

4750  LABEL  Temp$<0) 

4760  MOVE  125,4  !  DRAW  LABEL  2 

4770  LABEL  Temp$<l) 

4780  MOVE  130,4  !  DRAW  LABEL  3 

4790  LABEL  Temp$<2> 

4800  MOVE  130,53  !  LABEL  DIRECTIVITY  INDEX 

4810  LABEL  "  D  i  rec  t.  i  u  i  t  y  index  = " ;  DROUND  <  D  i  r  i  ndex ,  6  > 

4820  LDIR  0 

4830  !  SCALE  PLOTTING  AREA  TO  USER  DEFINED  UNITS 

4840  SCALE  -Dbf ul lscale,Dbful lscale,-Dbful lscale*.75,Dbful 1  seal e*. 75 

4850  !  DRAW  AXES  WITH  TICK  MARKS  EVERY  10  db 

4860  AXES  10,10,0,0,1,1 

4870  !  DRAW  DIRECTIVITY  PATTERN 

4880  FOR  Angl ecounter=0  TO  360-Resolution  STEP  Resolution 
4890  X=R ad i usual ues< Angl ecounter)*COS< Angle (Angl ecounter)+ 130) 

4900  Y  =  Radiusualues(Anglecounter-)*SIN(Angle(Anglecounter)  +  180) 

4910  IF  Ang 1 ecounter=0  THEN 

4920  MOVE  X , Y 

4930  ELSE 

4940  DRAW  X , Y 

4950  END  IF 

4960  NEXT  Anglecounter 

4970  !  COMPLETE  THE  DIRECTIVITY  PATTERN  PLOT  BY  DRAWING  A  LINE  FROM  THE 

4980  !  LAST  POINT  TO  THE  FIRST  POINT. 

4990  DRAW  Radi  usual ues < 0 ) *COS < Ang 1 e < 0 )  + 1 80 ) , Radi  usual ues<0)*SIN(Angl e(0)  +  130) 
5000  PEN  0 


5010  Yesno$="N" 

5020  INPUT  "MAKE  ANOTHER  PLOT  -  YxN.  DEFAULT  IS  <N)0",Yesno$ 

5030  IF  Yesno$=" Y"  THEN  4570 

5040  PLOTTER  7,5  IS  OFF 

5050  SUBEXIT 

5060  ELSE 

5070  SUBEXIT 

5080  END  IF 

5090  SUEEND 


TM  841152 


USING  THE  PROGRAM 

This  program  is  very  easy  to  use.  tost  input  responses  require  simple 
"yes"  or  "no"  answers.  The  process  of  digitizing  a  directivity  pattern 
consumes  the  most  time,  and  this  process  uses  the  stylus  of  the  graphics 
tablet  for  input.  To  execute  the  program,  the  operator  must  complete  the 
following  simple  steps: 

1.  Turn  on  the  HP  9845  computer 

2.  Turn  on  the  HP  9872A  plotter 

3.  Turn  on  the  HP  9111A  graphics  tablet 

4.  Load  "CALCDI"  from  the  proper  mass  storage  device 

5.  Press  the  "RUN"  key 

6.  Respond  accordingly  to  computer  generated  questions 
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Sample  wedge  pattern  generated  from 


previous  page 


Sample  spike  pattern  generated  from 


previous  page 


Sample  of  typical  pattern  generated 
from  previous  page  using  5°  resolution 


Sample  of  typical  pattern  generated. 


from  previous  page  using  2 


resolution 


